


> 

(N 



U. of MD PP #94-073 

DOE/ER/40762-023 

February 1, 2008 



Rapidity Gap in Jet Events at LEP200 



Hung Jung Lu 
t:!" ■ Department of Physics, University of Maryland 



College Park, Maryland 20742 



^vq I We analyze various perturbative mechanisms for the production of jet events 



containing rapidity gaps in e'''e~ annihilation at LEP200 energies. We found that the 
processes e~^e~ — > 7*7*, 7* Z, ZZ, WW -^ 4 jets generate gap events at an observable 
rate. We point out that LEP200 offers a unique opportunity for the study of gap 



o 

^r I events due to a smaller background of "fake-gaps" . 
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I. INTRODUCTION 

High-energy collisions involving hadronic final states typically create a large number of 
particles. A convenient kinematic quantity for the analysis of many-particle events is the 
rapidity variable 

where E is the energy of a final-state particle and P\\ its longitudinal momentum with respect 
to a reference axis. Rapidity is a natural longitudinal phase space variable. (See for instance 
IHJ^.) In fact, in e~^e~ annihilation into jets, one observes a roughly uniform distribution of 
hadrons per unit rapidity as measured with respect to the jet thrust axis, up to a maximum 
value ?/max ~ ln(-Ecm/"^7r)5 whcrc Ecra IS the total center-of-mass energy. A recent review on 
the use of the rapidity variable as a kinematic analysis tool can be found in Ref. 0J^ . 

In hadron-hadron collisions, the distribution of particles in the rapidity variable often 
contains gaps. That is, absence of particles in a rapidity interval. The existence of a gap 
may reflect the underlying exchange mechanism. For instance, a gap in diffractive scattering 
experiments can be interpreted as the exchange of a pomeron 0] . The possibility of using 
rapidity gaps as triggering signals in high-mass scale physics has motivated recent interest 
in the study of jet events containing gaps 0,^,0. 

In hadron-hadron collisions, the analysis of gap physics is complicated by the presence of 
spectator partons. These partons can interact with each other during the collision, creating 
additional particles which may contaminate and spoil a potential rapidity gap. The survival 
probability [0] of a rapidity gap is an important subject of study for unraveling of high-mass 
scale physics. 

In comparison, lepton-initiated collisions offer a cleaner environment, free of the uncer- 
tainties of survival probability considerations. The study of rapidity-gap events in lepton 
colliders can thus offer important complementary information to the study of gap events in 
hadron-hadron machines. 



Various perturbative mechanisms for generating events containing large rapidity gaps in 
e^e~ annihilation have been analyzed at low energies and at the Z-peak [§|,|^3 • Theoretical 
calculations indicate that these mechanisms generate gap events at an observable rate at 
the Z-peak. The perturbative QCD mechanisms are shown in Fig. |1]. Two color-singlet 
parton-pairs are produced in the final state. As these parton-pairs move away from each 
other, the color exchange between the two dipole systems is suppressed, and hadronization 
is dominantly confined to the phase space near each jet-pair. Hence we expect few or no 
particles to be produced in the rapidity interval between the two jet-pair systems. An 
interesting QED mechanism is shown in Fig. H, where a small-invariant mass photon decays 
into a high-rapidity hadron system, separated by a large rapidity gap from the other jet-pair. 
Although this QED mechanism is suppressed due to the smallness of Oem; its contribution can 
become important at low photon invariant-mass. Among all the QCD and QED mechanisms 
in Fig. |l] and Fig. H, the QCD qqqq mechanism of Fig. |l| (a) is the dominant channel 
for producing gap events at the Z-peak. As opposed to typical two-jet events in e+e~ 
annihilation which are distributed as 1 + cos^ 9 in the scattering angle of the thrust axis, the 
jet-pairs produced by the QCD mechanisms in Fig. [I| follow a distinctive sin^ 6 distribution. 

Random fiuctuations in the final hadron distribution can produce "fake gaps", which 
need to be distinguished from the gap events generated by the partonic mechanisms. The 
undergoing Monte Carlo studies are an important ingredient in the analysis of gap events. 

We encounter a unique experimental environment for studying gap events at LEP200, 
the future upgrade of the present e+e~ collider at CERN. The four-jet events from the decay 
of 7*7*, 7*Z, ZZ and WW boson pairs occur at a rate comparable to the two-jet events from 
of the decay of single 7* and Z bosons. Hence the "fake gap" background events at LEP200 
are expected to be reduced with respect to the situation at the Z-peak, where two- and 
three-jet events dominate over four-jet events. 

The study of gap events at LEP200 will also help to elucidate the nature of color fiow 
in quark fragmentation. For the decays of WW into four jets, Gustafson, Peterson and 
Zerwas [|lOl have proposed two possible scenarios. In the first scenario, the color fiow is 



confined to the quark-antiquark pairs of eacli gauge boson, and tlie two jet-pairs fragment 
independently. In the second scenario the string fragmentation occurs between the quark 
and antiquarks from opposite gauge bosons. A measurement of gap event cross section will 
therefore constitute an interesting test of the underlying fragmentation mechanism. 

In this paper we study the e~^e~ -^ 7*7*, 7*Z, ZZ, WW -^ 4 jets mechanisms for pro- 
ducing gap events at LEP200 energies. Due to the very different kinematic dependence of 
these channels (different flavor combinations, different invariant-mass resonant regions), we 
can neglect interference effects to first approximation, and consider the four mechanisms 
separately. (See also Ref. |ll[ for a discussion of interference effects in the production of 
heavy unstable particles.) The previously studied Z-peak mechanisms of Fig. |l| and Fig. |^ 
are suppressed by various kinematic factors (smaller coupling constant or large off-shellness 
in the gauge boson propagators), and do not have observable effects at LEP200. Our calcu- 
lations indicate that the dominant channels for producing gap events at LEP200 come from 
7*7* and 7*Z mechanisms. The WW mechanism has a smaller cross section, whereas the 
ZZ mechanism has a negligible contribution. 



II. 7*7* CASE 

The relevant diagrams are shown in Fig. ^. The final color-singlet jet-pairs [qiqi) (92^2) 
are effectively the decay products of the virtual photons. The dominant contribution of 
these amplitudes comes from the phase space regions where the photon propagators acquire 
large values, that is, near the photon poles. In fact, in the small photon invariant-mass 
limit, the cross section can be interpreted as the product of a{e'^e~ -^ 77) with the splitting 
functions of the two photons into respective qq pairs. Actually, color-singlet jet-pairs can 
also be formed by taking the quark from a virtual photon and the antiquark from the other 
virtual photon. That is, we can also form (^1^2) (^iQ'2) color-singlet pairs. However, for gap 
events the contribution from these "exchanged pairs" are strongly suppressed by the high 
virtualities carried by the photons, and also suppressed by an additional color factor 1/9. 
(This is analogous to the situation discussed in Ref. P] for QED mechanisms.) Hence, we 
need only to consider the contributions from the diagrams in Fig. |^. The expression for the 
differential cross section in the small photon invariant- mass limit is 

dxi a;^ + (1 — xi)^ dx2 X2 + (1 — X2Y 



dcos9 



tan2(0/2) + cot2(^/2) , (2.1) 



where Mi is the invariant mass of the (gigi) pair, M2 the invariant mass of the (^2^2) pair, 
Xi the longitudinal momentum fraction in the first jet-pair carried by gi, X2 the longitudinal 
momemtum fraction in the second jet-pair carried by q2, the scattering angle of the thrust 
axis, i/i the total center-of-mass energy, Qg^ and Qg^ the electric charge of gi and g2, 
«em = e'^/4:Tc, with e = gsinO^^ the electromagnetic coupling constant, g the weak coupling 
constant and ^w the weak angle, and finally, the triangular function is defined by 



A(s, M^, Ml) = ,Js^ + Mf + Mi - 2sMl - 2sM| - 2MlM^ . (2.2) 

The differential cross section is sharply peaked in the forward and backward beam direction. 
The apparent singularities at ^ ^ 0, tt in cot(^/2) and tan(^/2) can be traced back to 



the electron propagator in the t and u channeL These singularities are actually cut off 
by the photon invariant masses Mi and M2. For instance, for ^ ^ one can show that 
t -^ —56*^/4 — M1M2/S, hence an appropriate cut-off value of 6 is 

0,.. - '-^^. (2.3) 

s 

An analogous analysis for the u channel propagator leads to the upper limit ir — Ocut- There- 
fore the angular integral should be limited to the interval 9 G [Ocnt, tt — ^cut]- The final result 
is not very sensitive to the precise value of 6'cut because near the singularities the dependence 
of the angular integral on 6'cut is logarithmic: 

dar^f ^. (2.4) 



cut 



The cross section is also sharply peaked in the invariant mass variables Mf and M2. The 
apparent singularities in Mf, M2 — i> in Eq. (|2.1|) are physically cut off by the mass 
threshold of the qq systems. As in Ref. 0, we shall use the values M^, M^, M^, Mj,^ and 
My as representative values for the invariant mass thresholds for the uu, dd, ss, cc and bb 
systems. To obtain the total cross section, we integrate over all allowed flavor combinations, 
noting that allowed flavors depend on the integration variables M^ and M2. For higher 
values of M^ and M2 one has more allowed flavors. The upper bounds of Mf and M2 are 
given by the kinematic constraint A(s, Mf, M2) > in Eq. (|2.2| ) and also by the constraints 
requiring the existence of a rapidity gap. The integration limits for xi and X2 are set by the 
gap condition, too. (See Ref. |0,||.) 

In Fig. H we present the Lego plot of a typical 7*7* event. Notice that the two jet- 
pairs usually have a small invariant mass, generating thus a large rapidity gap in the central 
region. As in Ref. ||^ and Ref. [^ , we define g as the rapidity difference between the nearest 
jet centers from opposite jet-pairs. Since jet fragments typically are scattered over a radius 
~ 0.7 in the Lego plot 0, the effective gap is expected to be given by Qcs ~ (? — 1-4. Virtual 
photons can easily be transformed into vector mesons like p,uj,(j), J/ip and T. Hence, we also 
have the special cases where a vector meson recoils against a two-jet system or another vector 
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meson. In particular, the maximum value of rapidity gap is determined by events where 
we have two p mesons in the final state, that is, g = ln(s/M^). For a total center-of-mass 
energy oi E = ^/s = 180 [GeV], the maximum allowed gap is about 11 units of rapidity. In 
Fig. 1^ we present the integrated cross section of events having a value of nearest jet-center 
rapidity separation larger than g, for different energy levels. We have used aem ~ 1/128 
in our calculation. In Fig. ^ we plot the integrated cross section as function of energy, 
for different values of the gap-cut g. Given a total integrated luminosity of 500 [pb"^] per 



experiment expected for LEP200 ||T2[, the 7*7* mechanism should therefore generate gap 
events at an observable rate. 

Actually, the 7*7* mechanism also contributes at the Z-peak energy. In the small in- 
variant mass limit, there is no interference effects between this mechanism and the resonant 
QCD and QED mechanisms of Fig. |l] and Fig. ^ in the calculation of the total cross 
sections. (The interference terms vanish upon integration of the azimuthal angles 0i and 
02- See Appendix for the definition of (pi and 02- A similar cancellation was discussed in 
Ref. [^.) Therefore, each one of them can be obtained separately. In Fig. ^ we plot the 
integrated gap cross section for the QCD, QED and the 7*7* mechanisms at the Z peak 
{E = 91.17 [GeV].) We have used the results from Ref. 0,§], and imposed a representative 
invariant-mass cut of 30 [GeV] for the QCD and QED jet-pairs. We have used a Z-peak 
total visible cross section of 31 [nb]. The strong coupling constant is taken to be a^ ~ 0.13 in 
the calculations. As can be seen from the figure, the 7*7* mechanism has small contribution 
in the low-gap region, but becomes dominant in the large-gap region. Physically, the gap 
events generated by the 7*7* are distinct because the jet-pairs tend to be peaked towards 
the forward and backward directions, whereas in the cases of the resonant QCD and QED 
mechanisms, the scattering angle is distributed smoothly as sin^ 6 and 1 -|- cos^ 6. 



III. 7*Z CASE 

The relevant diagrams are shown in Fig. §. This mechanism is similar to the 7*7* 
case, but with a virtual photon replaced by a Z boson. As in the 7*7* case, we shall not 
consider (^152) and (^1^2) color singlet pairs, since the contribution of these "exchanged 
pairs" to gap events is suppressed by a large photon virtuality and also by the color factor. 
Due to the small invariant mass of the virtual photon and the large mass of the Z boson, 
we expect the decay fragments of 7* to be distributed in the large rapidity region on one 
side and those of Z in the lower rapidity region on the opposite side. In general we will 
encounter a lopsided gap as depicted in Fig. |^. A special case would be a vector meson 
recoiling against the two-jet system from the Z boson. The maximum value of rapidity 
gap g = ln[(s — M'^)/MpMz] is determined by events where we have a p meson recoiling 
against two jets that are distributed symmetrically with respect to the thrust axis. For a 
total center-of-mass energy of E = ^/s = 180 [GeV], the maximum allowed value of g is 
about 5.8 units of rapidity. 

The differential cross section in the limit of small photon invariant mass is given by 

'^''~ 877 ^"^ Ml {Ml - Miy + {TzMzYis - M^y^^''^^^''^^^'' 
dxi dx2 d cos el (1 + Ma/s^) [xl + (1 - xif 

' {QfQl' + Qe 'Q"/) {4 oot'{9/2) + (1 - x,f tan2(^/2) 
+ (gf Q^/ + Q]:'qI') {xlt^n\9/2) + (1 - a:2)^cot2(^/2) 



M 



2 



+4q:q^,^ [xi + (1 - x,y\ X2(i - X2)y (3.1) 

In the above expression, ttem = e^/47r = g^ sin^ 6'w/47r is the electromagnetic coupling 
constant and a-w = 5'^/4vr the weak coupling constant. Mi is the invariant mass of the 
(QiQi) pair, M2 the invariant mass of the (q'2^2) pair, Mz and Tz the mass and width of 
Z boson, Xl the longitudinal momentum fraction of the virtual photon carried by gi and 
X2 the longitudinal momentum fraction of the Z boson carried by q2 (see Appendix), 6 the 
scattering angle of the thrust axis, ^/s the total center-of-mass energy, Qg-^ the electric charge 



of qi, Qf, Qg , Q^ and Q^^ the weak charges of e and q2, where for a fermion of isospin // 
and electric charge Qf we have 



n — 1*^/ 1 — [ sec^w-^/ — sin6'wtan6'iy<5/\ 
^ \Qf) " V - sin ^w tan O^Qf ) ' 

Q'f = Qf + Qf. (3.2) 



Due to the narrow width of the Z boson, we can make the approximation 



/ 



dMi TT 



(3.3) 



{Ml - M^z? + iXzMzY VzMz 
and replace M"^ by M"^ elsewhere in Eq. ( p.l| ). 

The 9 angle dependence in Eq. ( ^.11 ) implies the existence of a forward-backward asymme- 
try in the cross section, which is also clear from the explicit dependence on the left-handed 
and right-handed weak charges. The singularities at ^ 0,7r in cot(0/2) and tan(^/2) 
can be traced back to the electron propagator in the t and u channel. These singularities 
are physically cut off by the photon invariant mass Mi. For 0^0 one can show that 
t ^ — (s — M|)^^/4 — MfM|/(s — M|), therefore an appropriate cut-off angle is 

Similarly, by analyzing the u variable, one can show that the upper limit for the ft integration 
should be TT — 6'cut- Define 

Ig{M^)= r "'rfcos0tan2(0/2)= r "' t/cos0cot2(0/2), (3.5) 

the expression of the cross section in Eq. ( p.l| ) simplifies to 



"e-Cq.^q, 1^(5 -M|)2 

^A(.,M,^M|) 

dxi x^ -|- (1 — xiY 

rfx2 {(l + M'z/s') [xl + (1 - x^f] Ie{Ml) 

+ 8(M|/s)x2(l-X2)}. (3.6) 
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To obtain the total cross section, we sum over all allowed flavor combinations. As in the 
7*7* case, we use the values M'^^M'^^M'^^M'^ij^ and M^ as representative values for the 
invariant mass thresholds for the uu,dd,ss,cc and bb systems. The xi, X2 integrations, as 
well as the Mf integration, are also limited by the gap condition 0,§|. In Fig. [1^ we plot 
the integrated cross section of events having a value of nearest jet-center rapidity separation 
larger than g, for different energy levels. We have used aw ~ 1/29.8 and sin^ 6'w = 0.233 



in our calculation. In Fig. |TI| we plot the integrated cross section as function of energy, for 
different values of the gap-cut g. For the projected integrated luminosity of 500 [pb~^] per 
experiment at LEP200, the '~f*Z mechanism therefore also generates an observable number 
of gap events. For smaller values oi g {g < 4) the 7*Z mechanism dominates over the 7*7* 
mechanism. As discussed previously, the 7*Z mechanism cannot not generate gap events 
with g larger than ~ 5.8, hence for large gaps {g > 5) the 7*7* mechanism dominates. 
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IV. ZZ CASE 

As opposed to the 7*7* and ^*Z mechanisms, the direct jet-pairs from the decay of ZZ 
and WW particles are not capable of generating events with large rapidity gaps. This is 
because at LEP200 the ZZ and WW pairs are produced right above the threshold, thus 
the two gauge bosons recoil slowly against each other. For instance, at a center-of-mass 
energy of 200 [GeV], we have only 1.4 units of rapidity separating two recoiling W particles, 
which means in practice no observable rapidity gap. However, the "exchanged pairs" as 
those shown in Fig. |T^ provide a feasible channel of producing gap events. It is interesting 
to recall that the exchanged-pair contributions are negligible in the 7*7* and 7*Z cases due 
to the large invariant mass carried by the virtual photons, therefore the direct pairs are 
responsible for gap events there. The situation is exactly reversed for the ZZ and WW 
cases, where the exchanged pairs are the responsible mechanism. 

Let us consider first the case of a ZZ pair produced exactly at threshold. In the center- 
of-mass frame, each of these static Z particles decays into a back-to-back quark-antiquark 
pair. In the phase space where q\ becomes coUinear with g2 (hence q\ becomes also collinear 
with ^2), we can form small-invariant-mass color-singlet jet pairs and have a rapidity gap. 
Naturally, to have an observable cross section, the center-of-mass energy should be larger 
than the threshold value. As the center-of-mass energy becomes higher and higher, the 
crossing of a quark from the phase space of a Z boson to the phase space of the opposite Z 
boson becomes less and less likely. Therefore, we expect the gap event cross section to rise 
from zero to a maximum value somewhere above the threshold, and then to decrease with 
increasing center-of-mass energy. 

The crossing of particles in the kinematic phase space substantially complicates the 
calculation of the gap event cross section. This is because the gap condition becomes a 
unnatural cut in the integration of the various kinematic variables. In particular, as opposed 
to the 7*7* case and the 7*Z case, the azimuthal angles 0i and 02 (see Appendix) can no 
longer be integrated out analytically. 
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We compute the helicity amplitudes in terms of spinor products defined in Ref. [|T3 . 
After using the narrow width approximation for the Z boson propagators to integrated out 
the jet-pair invariant masses Ml and M2^ we obtain 



„4 



^^ = „^ooT^ ?# — ^ A(s, M|,M|) dxidx2d(t)id(t)2dcose 



7r(32r7Mz s) 

E E Qf''QfyQf,''^"\nHe,H,,,H,,)\\ (4.1) 

where the first sum is over all allowed quark flavors, and the second sum is over all helicities 
of the electron and the final quarks gi and g2- We have included in the formula an 1/2 
factor to account for the double counting of flavors in the summation of gi and g2- The 
various charges in Eq. ( |4.1j ) are the weak charges of the electron, gi and q2, as defined in the 
previous section. (We use indifferently {R,L) and (+, — ) to denote right and left helicity.) 
There are eight helicity amplitudes, but half of them are needed since the other half can be 
obtained by simple helicity conjugation: 

Ti+ + +) = J {kiqi]{k2q2) {-[kiq2]{kiqi) + [qiq2]{qiqi)} 

4 _ ( _ _ ^ 

H — [kiq2]{k2qi) \-[kiqi]{kiq2) + [qiq2]{q2q2) \, 

4 r _ _ . 

T{+ + -) = - [kiqi]{k2q2) [-[kiq2]{kiqi) + [gig2](gigi)| 

[kiq2]{k2qi) \ [kiqi]{kiq2) + [qiq2]{q2q2) \ , 

T{+ -+) = -- [kiqi]{k2q2) {[kiq2]{kiqi) + [gig2](gigi)| 

4 r _ _ _ _ -, 

+ - [kiq2]{k2qi) [-[kiqi]{kiq2) + [qiq2]{q2q2) \ , 

T{+ --) = J [kiqi]{k2q2) {-[kiq2]{kiqi) + [q2qi]{qiqi) \ 

4 r -1 

+ - [kiq2]{k2qi) {-[kiqi]{kiq2) + [q2qi]{q2q2)\- (4.2) 

The spinor products in these amplitudes can be evaluated numerically with the formulas 
given in the Appendix. 

Due to the crossing of quark pairs, the thrust axis in general is not aligned with the 
gauge boson direction. We will denote the scattering angle of the thrust axis by O, which 
in general differs from the gauge boson 6 angle that appears in the expression of the cross 
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section of Eq. (|4.1|) . In the Monte Carlo integration of the cross section in Eq. |4.1| , we only 
retain those events containing a rapidity gap larger than the gap cut g. The result is plotted 



in Fig. |T3|. We notice that large gap events from the ZZ mechanism have cross sections 
much smaller than 0.01 [pb]. Hence, this mechanism is not expected to have observable 
effects at the projected luminosity of LEP200. In Fig. IT^ we plot the angular distribution 
of the cross section for E = 192 [GeV] and g = 3. We see that the cross section has two 
smooth humps, as opposed to the 7*7* and 7*Z cases where the cross section has two sharp 
peaks ~ 1/6 in the forward and backward beam directions. 



V. WW CASE 



The kinematics of the WW case is similar to the ZZ case. Diagrammatically, the 
production of WW boson pairs is interesting because it involves the triple-gauge-boson 
vertex of the electroweak theory. The Feynman diagrams of this mechanism is shown in 
Fig. |T5|. The Ui, U2, and di, d2 quarks represents any up-type (m, c) and down-type {d,s,b) 
quarks and antiquarks. As explained in the previous section, at LEP200 the direct-pair 
contributions are not capable of producing events with large rapidity gaps, therefore we 
need only to focus on the exchanged-pair contributions of Fig. |15[ Like in the ZZ case, we 
also expect the gap event cross section to reach a maximum somewhere above the threshold, 
since at higher energies there is less phase space for the formation of exchanged pairs. 

The differential cross section for the WW case is given by 



da 



a 



w 



5127i{TzMzsy 



A(s,M^,M^) dxidx2d(j)id(j)2dcos9 



T^ 



s s- Ml + iVzMz 



sin^ 9- 



w 



+ 



'sin^^ 



w 



+ 



i - sin^ e- 



w 



T{ 



-To 



(5.1) 



s-Ml + iVzMz) ' 2i 

where M^i and Tz are the mass and width of the W boson, A as defined in Eq. ( ^.21) . 
We have summed over all allowed quark flavors and made following approximation for the 
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summation over the Cabbibo-Kobayashi-Maskawa matrix elements (See for instance Ref. 



0) 



E \v., 



2. 



(5.2) 



j^d,s,b 



The T^ and T^ amphtudes come from the s-channel diagram in Fig. |1^ (a), and are given 
by 

T+ = 4 [diU2]{d2Ui) \^[kiUi]{uik2) + [kidi]{dik2) 
+ 4 [kiU2]{d2k2) [[diU2]{u2Ui) + [did2]{d2Ui) 

- 4 [kidi\{uik2) [[u2Ui\{uid2) + [u2di]{did2) 
T~ = 4 [diU2]{d2Ui) [[k2Ui]{uiki) + [k2di]{diki) 

+ 4 [k2U2]{d2ki) \^[diU2]{u2Ui) + [did2]{d2Ui) 

- 4 [k2di\{uiki) \[u2Ui]{uid2) + [u2di]{did2)j, (5.3) 

where we have used the symbol ui, di, U2 and ^2 to denote the final state quarks as well as 
their respective final momenta. The T^ amplitude is given by 



T2 = 4 [uiki]{k2U2) {[uidi]{d2Ui) - [kidi]{d2ki)y 



(5.4) 



The spinor products in these expressions can be evaluated with the formulas given in the 
Appendix. 

As in the ZZ case, the thrust axis angle 6 in general differs from the gauge boson angle 
6 in the expression of the cross section of Eq. (|5.1| ). In the Monte Carlo integration of the 



cross section in Eq. ( ^.l]) we only retain those events with a rapidity gap larger than the 
gap cut value g. The resulting cross section is given in Fig. |16[ We see from the plot that 
the WW mechanism can contribute somewhat to the production of gap events at LEP200. 
For gap cuts g between 3 and 4 units of rapidity, the maximum cross section is reached 
around 5 to 10 [GeV] above the threshold. In Fig. |1^ we plot the angular distribution of 
the cross section for E = 170 [GeV] and g = 3. As in the ZZ case, the cross section shows a 
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double-hump structure, with moderate depletion near the beam direction and the direction 
perpendicular to it. 

VI. CONCLUSION 

We have analyzed in this paper various perturbative mechanisms for the decay of elec- 
troweak gauge bosons into jet events containing a rapidity gap. At LEP200 energies the 
dominant contributions come from the 7*7* and 7*Z mechanisms. To fix an idea, for a gap 
cut g = 4:, '~f*Z is the dominant contribution with a cross section a ~ 0.04 [pb], and the 7*7* 
mechanism contributes another a ~ 0.015 [pb]. For smaller gap cuts the cross section is 
larger. However, if the gap cut is too small, the background "fake gap" events from random 
fiuctuation of final-state hadron fragments can be significant. For the 7*7* events the jet 
fragment distribution typically is peaked towards the forward and backward beam direction. 
For the '~f*Z process one would typically observe a forward or backward jet system recoiling 
against a low-rapidity two-jet system on the opposite side, hence have a lopsided rapidity 
gap. It is also likely to observe the special cases where the virtual photons are converted 
into a vector mesons instead of jet-pairs. 

The WW mechanism may also contribute to the production of gap events. Our analysis 
indicates that the maximum event cross section is reached 5 to 10 [GeV] above the threshold. 
These events have a relatively smooth distribution in the scattering angle. The measurement 
the experimental gap event cross section here will help to test the hypothesis that the color- 
fiow is confined to the collinear jets. Finally, the ZZ mechanism has a very small cross 
section at LEP200, and we do not expect it to contribute in practice. 

We thank Stanley J. Brodsky, Valery A. Khoze and Joseph Milana for helpful conversa- 
tions. This work is supported by Department of Energy contract DE-FG02-93ER-40762. 
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APPENDIX A: FOUR-JET KINEMATICS 

All the production mechanisms in this paper involve the decay of two gauge bosons into 
four final-state quarks. The graphic representation of the particle momenta is shown in Fig. 
|TH|. We denote the momenta of the incoming positron and electron by ki and k2, respectively, 
and we will use qi,qi,q2, ^2 to denote the final-state quarks as well as their momenta. The 
momenta carried by the first and second gauge bosons are denoted by Pi and P2, Therefore 
we have 



-Pi = qi + qi, 

P2 = q2 + q2- 



(Al) 



That is, the (gigi) pair is the decay products of the first gauge boson, and the (52^2) pair is 
the decay products of the second gauge boson. 

At the center-of-mass frame, the momenta of the four quark jets can be parametrized in 
terms of seven kinematic variables: M^, M|,xi,X2, ^, 0i, 02? where Mf = P^,M2 = P2 are 
the invariant masses of the (qiqi) and (q'2^2) pairs, Xi and X2 the longitudinal momentum 
fractions of gi within the (gigi) pair and q2 within the (52^2), (^ is the gauge-boson scat- 
tering angle, and 0i,02 the azimuthal angles of gi and ^2 with respect to the gauge boson 
direction. Notice that 6 coincides with the thrust axis scattering angle for the 7*7* and 'y*Z 
mechanisms. However, for the ZZ and WW mechanisms the thrust angle B defined by the 
direction of gi + g2 is different from 6 due to the crossing of jet-pairs. The following is the 
explicit representation of the initial and final four-vectors: 



ki 
k2 

<12 



2 

VI 

2 
1 

2 
1 



1; k 



1; -k 



xiP^ + xiP^] 2v/^^Mi pi^ + (X1P1+ - XiPf ) P 
xiP^ + xiP^; -2^xiXiMi pi^ + (xiPi+ - xiPf ) P 



X2P2+ + X2P2"; 2y/^^M2 P2± - (a;2P2+ - X2P2 ) P 
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g2 = 2 [ ^2 A^ + X2P2 ; -2^/x^M2 P2± - {,X2P2 - X2P2 



(A2) 



It is convenient to choose the z axis along the direction of the momentum of the first gauge 
boson Pi, This will facilitate the calculation of spinor products later. Using the coordinate 



system shown in Fig. |18|, we have 



k = cos 6z — sin 0x, 
Pi_L = cos 0ix + sin 0iy, 
P2± = cos 02X + sin (j)2y, 

P = z. 

The quantities Xi, Xi, X2, a;2, P^, Pi, P2 , P2 satisfy the following constraints 

Xi +Xi = 1, 

X2 + X2 = 1, 

P^P{ = Ml 
P2+P2- = Ml 



(A3) 



(A4) 



From conservation of energy and momentum, we can obtain Pi, P^ , Pi, P2 explicitly in 
terms of Mf , M2 and s 

1 



^1 
P 



2Vi 
1 



2Vi 
1 



' "2v/i 
1 



2Vi 



s + M^- Ml + A(s, Ml Ml) 
s + Ml- Ml - A(s, Ml Ml) 
S-MI + Ml + A(s, Ml Ml) 
s-Ml + Ml - A(s, Ml Ml) 



(A5) 



where A(s, Ml Ml) is the triangular function defined in Eq. ( p.2| ). 

The spinor products (see Ref. [|13| for a review on the helicity method) between the 
various four-vectors can be computed explicitly. We have 
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{kih) 
{hqi) 
{hqi) 
{hq2) 
{hq2) 
{hqi) 
{hqi) 
{k2q2) 
{k2q2) 

{qiqi) 
(9292) 
{qiq2) 
{qiq2) 
{qiq2) 
{qih) 



-vs, 

\^P2X2 
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\^P^Xi 
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X1X2P1P2 
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2 

1/2 e 

sin - - 
2 

1/2 e 

sin - - 
2 

1/2 ^ 
sin - - 
2 

1/2 ^ 
cos - — 
2 
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2 



v^PfaJi 

\^P^X2 



1/2 ^ 



COS - e 
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cos-e*^S 

1/2 9 • . 
cos- e*^^ 
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«<P2 



*<P2 



1/2 



g*0l _ 



1/2 
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1/2 



1/2 



gi<^i 



e*0i _ 



a:iX2Pi+P2+ 

X,X2P^P+ 
XiX2P^P^ 
XiX2P^P^ 



1/2 



:,i</'2 



1/2 
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1/2 



1/2 



3«02 



(A6) 



All other spinor products can be obtained by the antisymmetry operation or by helicity 
conjugation: 



(ab) = —{ha) 
[ah] = -[ha] = {hay 



(A7) 
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FIGURES 
FIG. 1. Perturbative QCD mechanisms for generating rapidity-gap events at the Z peak. The 

dashed hnes indicate that the produced partons are in color-singlet state, (a) Two final-state 

quark-antiquark pairs, (b) A quark-antiquark jet pair and a two-gluon jet pair. 

FIG. 2. A QED mechanism for generating rapidity-gap events at the Z peak. The dashed lines 
indicate that the produced partons are in color-singlet state. 

FIG. 3. The e+e^ -^ 1*1* ~^ QQQQ mechanism for generating gap events. 

FIG. 4. A typical Lego plot for a gap event from the 7*7* mechanism. The physically observed 
gap is expected to be ^eff — 9 — 1-4 due to the spreading of hadron fragments around each jet. 

FIG. 5. Gap event cross section for the 7*7* mechanism, as function of the gap-cut g and for 
two different values of energies. 

FIG. 6. Gap event cross section for the 7*7* mechanism, as function of the energy and for 
different values of the gap-cut g. 

FIG. 7. The 7*7* contribution to the gap-event cross section at the Z-peak (solid line). Also 
plotted are the contributions from the QCD (dashed line) and QED (dotted line) mechanisms (see 
Fig. I, Fig. I and also Ref. g.) 

FIG. 8. The e"'"e~ — > 7*Z -^ qqqq mechanism for generating gap events at LEP200. 

FIG. 9. A typical Lego plot for a gap event from the 7*Z mechanism. Notice that the rapidity 
gap is shifted towards the virtual photon side. 

FIG. 10. Gap event cross section for the 7*Z mechanism, as function of the gap-cut g and for 
two different values of energies. 
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FIG. 11. Gap event cross section for the 7*Z mechanism, as function of the energy and for 
different values of the gap-cut g. 

FIG. 12. The e~^e~ -^ ZZ -^ QQQQ mechanism for generating gap events at LEP200. Notice 
the exchange of quarks and antiquarks from different gauge bosons to form color-singlet jet-pairs 

FIG. 13. Gap event cross section for the ZZ mechanism, as function of the energy and for 
different values of the gap-cut g. 

FIG. 14. Angular distribution of gap event cross section for the ZZ mechanism, where is the 
thrust axis scattering angle. 

FIG. 15. The e^e~ -^ WW -^ qqqq mechanism for generating gap events at LEP200. Notice 
the exchange of quarks and antiquarks from different gauge bosons to form color-singlet jet-pairs 

FIG. 16. Gap event cross section for the WW mechanism, as function of the energy and for 
different values of the gap-cut g. 

FIG. 17. Angular distribution of gap event cross section for the WW mechanism, where Q is 
the thrust axis scattering angle. 

FIG. 18. Initial and final momenta involved in e^e~ -^ qqqq- 
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